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Abstract 

Quantum Hamilton-Jacobi Theory and supersymmetric quantum mechanics (SUSYQM) are two par- 
allel methods to determine the spectra of a quantum mechanical systems without solving the Schrodinger 
' equation. It was recently shown that the shape invariance, which is an integrability condition in SUSYQM 

O 

' formalism, can be utilized to develop an iterative algorithm to determine the quantum momentum func- 

' tions. In this paper, we show that shape invariance also suffices to determine the eigenvalues in Quantum 

^ ' Hamilton-Jacobi Theory. 

p<j ' Quantum Hamilton-Jacobi Theory and supersymmetric quantum mechanics (SUSYQM) are two very 

different methods that give eigenvalues for quantum mechanical systems without solving the Schrodinger 
differential eigenvalue equation. 



Supersymmetric quantum mechanics is a generalization of Dirac's ladder operator methoqj for the har- 
,■ monic oscillator. This method consists of factorizing Schrodinger's second order differential operator into 



two first order differential operators that play roles analogous to ladder operators. If the interaction of a 



\ quantum mechanical system is described by shape invariant potentials [HH], SUSYQM allows one to generate 

lO i all eigenvalues and eigenfunctions through algebraic methods. 

^SJ ' Another formulation of quantum mechanics, the Quantum Hamilton-Jacobi (QHJ) formalism, was devel- 

OO ' oped by Leacock and Padgett [3] and independently by Gozzi 4 . It was made popular by a series of papers 

o ; 

I by Kapoor et. al. [5]. In this formalism one works with the quantum momentum frmction (QMF) p{x), 
which is related to the wave function ip through the relationship p{x) — —ipf{x)/ilj{x), where prime denotes 
differentiation with respect to x. Our definition of QMF's differs by a factor of i = from that of ref. 

^ ' [3l[5l[6], where they define p{x) = —i ^j^^^ ; we use p{x) — — ^^j^^^ . It was shown, on a case by case basis, 
that the singularity structure of the function p{x) determines the eigenvalues of the Hamiltonian [31 131 [5] 
for all known solvable potentials. Kapoor and his collaborators have shown that the QHJ formalism can be 
used not only to determine the eigenvalues of the Hamiltonian of the system, but also its eigenfunctions [7] . 
They have also used QHJ to analyze Quasi-exactly solvable systems where only an incomplete set of the 
eigenspectra can be derived analytically and also to study periodic potentials [8]. It is important to note 
that all cases worked out in Refs. [31 E] satisfied the integrability condition known as the translational shape 
invariance for which the SUSYQM method always gave the exact result. 



'^e-mail; ccherqu@luc.edu 
■^e-mail: yev@freeshell.org 
^e-mail: agangop@luc.edu 
■^Dirac attributes this method to Fock. 
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In an earlier paper [B], one of the authors (AG) and his collaborators had connected the Quantum 
Hamilton- Jacobi formalism to supersymmetry (SUSY) and shape invariance. Using shape invariance, it was 
shown that QMFs corresponding to different energies are connected via fractional linear transformations, 
and a general recursion formula was given for quantum momenta of any energy. Since QMFs are related 
to wave functions, Ref. provided an alternate recursive procedure to derive eigenfunction directly from 
shape invariance. 

In this paper, we show that shape invariance also suffices to determine the singularity structure of the 
QMFs, and hence provide an explanation why translational shape invariance also determines the eigenvalues 
of the system with the QHJ formalism. 

In the QHJ formalism the spectrum of a quantum mechanical system is determined by the solution of 
the equation: 

p'^{x,ao) - p'{x,aa) ^V{x,ao) - E (1) 

Here we have chosen our units such that h — 1 and 2m =1. ao is a parameter characterizing the strength 
of the potential. This equation is related to the Schrodinger equation 

+ {V{x,ao)-E)i^ = (2) 

via the correspondence 

P= — (-^] whence tlj{x) ^ e~ -f ^^^^"^^ . (3) 
\V J 

In SUSYQM, the supersymmetric partner potentials V±{x,ao) are given by W'^{x,ao) ± W'{x,aQ) respec- 
tively. VF(x, ao) is a real function, called the superpotential [Ij. In this paper V-{x,aQ) is chosen as the 
potential for a Hamiltonian H^(x, ao) whose eigenvalues we will be seeking. Thus, in terms of the superpo- 
tential W{x, ao), the Quantum Hamilton- Jacobi equation we will be solving is 

p^{x,ao) - p'{x,ao) = W^{x,ao) - W'{x,ao) - E . (4) 

We will assume that our superpotential is such that SUSY remains unbroken for a range of values of 
the parameter ao. This implies that the ground state eigenvalue E^ ^ = [1], and the correspond- 
ing eigenfunction ijj'^ {x) ~ g- / ^(^'"o)da; normalizable. By differentiating ip^^ one can show that 
VF(a;, ao) = — V'o ''(2^) /^o \^)- ^ comparison with Eq. [3] shows that we should expect the following 
limit: 

p [x) = W{x, ao) for E ^ 0. (5) 

(Note that the partner Hamiltonian i/_|_(a;, ao), built using the potential V+{x, ao), has the same set of 
eigenvalues Er^^ — ^L+u except for the groundstate.) 

As stated earlier, in this paper our focus will be on i7_(x, ao). In SUSYQM literature, it is customary 
to denote the eigenfunctions of H- by and the corresponding QMF by p-{x) = — t/;'^' y?/;'^). Since 
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we are going to be working with and its solutions exclusively, we will henceforth suppress the "-" 
subscript /superscript . 

The QHJ method is centered around analyzing the pole structure of the quantum momentum function 
p{x). Since the n-th eigenfunction of the Hamiltonian has n nodes within the domain of the potential, the 
corresponding momentum function p{x) will have n singular points which can be shown to be simple poles 
[5], each with a residue of —1. These are known as the "moving poles". Hence, a contour integral of p{x) 
around the moving poles of the system yields: 

p{x) dx = —2TTin . (6) 

This is the quantization condition for QHJ. For the sake of clarity, let us reiterate that the factor of — i in 
the above quantization condition is due to our definition of the QMF differing from that of refs. O [5l [6] 
by a factor of i. Our choice is guided by our desire to keep p{x) close in its definition to the superpotential 
W{x) so Eq. ^ holds, without any factor of i. For n ~ 0, p{x) would be identical to W{x). 

The strength of QHJ derives from the fact that the contour of the above integration can be deformed 
to enclose "fixed poles" instead (albeit traversing in the opposite direction [3 [5]). Thus, the quantization 
condition can now be written as 



p{x) dx 



-27rm . (7) 



p{x) dx 

origin 

The difference in sign between Eqs. ([6]) and ^ is due to the change in direction in which the integration is 
carried out. In this paper, we only consider systems that have fixed poles at the origin and/or oo. 

We now consider the impact of shape invariance on the Hamiltonian- Jacobi formalism. Shape invariance 
implies the following relation among partner potentials of supersymmetric quantum mechanics: 

V+(r,ao) = V-{r,ai)+R{ao), 

i.e., potentials V+(r) and V-{r) have the same functional dependence on the coordinate variable r and only 
differ from each other in the value of the parameter ao and an additive constant R{ao). For all known shape 
invariant potentials, it is possible to express the constant R{ao) as a difference: 

R{ao) = 9{ai) - g{ao) ■ (8) 

We assume that i?(ao) can be written as a difference of ^-functions for all shape invariant problems. 
The shape invariance condition, when written in terms of the superpotential W, takes the form: 

T,r9/ X dW(r,ao) ^,^9, , dW(r,ai) , . . 

dr dr 

We will consider only the additive shape invariance with ai = ao + 1, and use it to determine the structure 
of the superpotential near zero and 00, two assumed end points of the domain. It is essential that the 
knowledge of the superpotential at these boundary points be known in order to determine eigenvalues using 
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QHJ. In particular, we will not consider cyclic shape invariant systems [51 llOj because their superpotentials 
have an oscillatory behavior at large distances; i.e, they are not well defined at infinity. 

Here we shall divide all superpotentials into two categories: those with infinite domain (0, oo) or (—00, 00) 
and those with finite domain, such as (xi,X2). Superpotentials that are given in algebraic form or those 
in which a transformation such as r = converts it into an algebraic form will be considered in this 
paper. These potentials are defined over an infinite domain. We call them as "Algebraic" and "Hyperbolic" 
respectively. In the following, we will show how the shape invariance suffices to determine eigenvalues of 
generic algebraic and hyperbolic potentials. (We will not be considering superpotentials with finite domain 
which we define as trigonometric.) 

Algebraic Potentials: 

In this category, we assume that the potential is explicitly expressed in an algebraic form. Let us assume 
the superpotential has the following structure near origin: 

W(r,ao) = bkiao)r'' = ^^1^ + bo{ao)r° + h{ao)r^ + b2{ao)r^ + ••• 

k=-l 

We have not included a term more singular than 1/r in the superpotential to avoid the particle's "falling 
to the center." Particles cannot be sustained in a bound state in an attractive singular potential if the 
singularity is higher than —472 [H]- Additionally, we are interested in the simple poles of the QMF at 
boundary points, which are expected to be related to the same order singularities of the superpotential at 
those points. 

The derivative of the superpotential is then given by 

dW{r,ao) &-i(ao) , n , h t \ , / \ , 

= 5 ^ + bi[ao) + 2b2{ao)r + ■ ■ ■ ; 

ar 

and the square of the superpotential is given by 

W^{r,a^) = (6_i(ao))'r-2 + (2 6_i(ao) 6o(ao)) + ((6o(ao))' + 2 6-i(ao) &i(«o)) + ••• 

Substituting these expansions into the shape invariance condition ^ , we arrive at the following constraints 
on coefficients of the superpotential W{x, a^): 

from the r^^ terms (6_i(ao))^ ~ 6-i(ao) ~ (^-1(0:1))^ + ^-i(«i) 

from the r^^ terms (2 6_i(q!o) &o(ao)) = (2&-i(q:i) 6o(«i)) 

/ 2 \ (10) 

from the r° terms ( (6o(ao)) +26_i(ao) ^i(ao)) +&i(ao) 

= ((5o(ai))' + 26_i(ai) 6i(ai)) - 6i(ai) + i?(ao) 

The first of these three constraints of Eq. PT7)) is a quadratic equation for the coefficient 6_i(ai) in terms 
of 6_i(ao)- It has two solutions: 6_i(ai) = — 6_i(ao) or 6_i(ai) — 6_i(ao) — 1. Since 6_i(ao) is the 
coefficient of the most dominant term and hence determines the structure near the origin, its value plays an 
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essential role in deciding whether SUSY is unbrokerO. Hence, the first solution 6_i(ai) = — 6_i(q!o), is not 
acceptable if both 6_i(q!o) and 6„i(q!i) were to keep the system in a parameter domain needed for unbroken 
SUSY. Thus, the first constraint (jTU]) implies the relationship 6_i(q;i) — 6_i(q;o) — 1- For ai = ao + 1, this 
gives the difference equation 6_i(ao + 1) = 6_i(q;o) — 1, whose solution is given by 

6_i(q;o) = — ao + constant; (11) 

we choose the constant = 0. This choice, while it simplifies the calculation, has no effect on the final result. 
The second constraint leads to the difference equation 

ao boiao) = ai 6o(ai) . (12) 

In other words, the product ao 6o(c«o) does not depend on the parameter ao- Hence, we have 

boiao) = — , 

where /? is a constant independent of the parameter ag. Thus, near the origin, the structure of the superpo- 
tential Vl^(a;,ao) is given by: 

W(-.«o)Ui„ = -V + ^ + --- (13) 

Now, let us explore the structure of the superpotential near oo. To do this, we define a variable u = ^, and 

W {l/u , ao) = W (u , ao). The shape invariance condition of Eq. ((9]) thus transforms to 

^r? / N y dW (u, ao) ~5 , , ^ dW (u , ai) . 

[u , ao) - \ ' =W^{u, ai) + \ ' + i?(ao). (14) 

au au 

Analogous to the expansion near the origin, let us expand W (i , ao) in powers of u as we have done earlier. 

°° Q i(ao) 

W{u, ao) = Cfc(ao)u'' = — — — + co{ao)u° + ci{ao)u^ + C2(ao)M^ + ••• 

u 

k=-l 

w2/ ^ (c_i(ao))^ (2c_i(ao) co(ao)) f , ^^2 . „ , ^ / ^\ , 
W (m, ao) = 5 1 h (co(ao)) + 2 c_i(ao) ci(ao) H 

dW(r,ao) n dW(u,ao) , . , ■? „ / n 

J^ = -" = c_i(ao) - ci(ao)w2 _ 2c2(ao)u3 + ... ; 

dr au 

Substituting these expressions in the transformed shape invariance condition of Eq. (|14p . we get 

from the terms (c_i(ao))^ = (c_i(ai))^ ; (15) 

from the r^-'^ terms (2 c_i(ao) co(ao)) (2 c_i(ai) co(ai)) ; and (16) 

from the r" terms (co(ao))^ + 2c_i(ao) ci(ao) +C-i(ao) 

= (co(ai))^ + 2c_i(ai) Ci(ai) - c_i(ai) + i?(ao) . (17) 



^The limiting values of the superpotential at two boundary points of the domain determine whether SUSY is broken or 
unbroken. If the value of the superpotential at both end points have the same sign SUSY is broken and SUSY is unbroken if 
both end points have opposite signs. Thus, if the sign of W at one end is altered relative to the other end, the SUSY goes from 
being unbroken to broken or vice versa. For more details, see Ref. [Tl|2] 
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From Eq. psp . we get c_i(q;i) = ±c_i(ao). As stated earlier we assumed that SUSY remains unbroken 
for a range of values of the parameter, and since c_i(q;i) — — c_i(q;o) would lead to breaking of SUSY, 
we choose c^i{ai) = c_i(ao) = c_i; i.e., this coefficient does not depend on the parameter ctQ. This 
then further leads to co(q!o) — co{ai) — cq, another constant. Eq. (|17p now yields, 2c_i ci(q!o) + c_i = 
2 c_i ci{ai) — c-i + g{ai) — g{ao). Which gives 



We can write this equation as 



, N ,5(01) / ^ , g("o) , 1 /1o^ 

2c_i 2c_i 



/(ai)=/(ao) + l , (19) 



where /(ao) = ci(q;o) + • Since ai ~ ao + 1, the solution of Eq. is /(ao) = ao + A. Hence, we 

have ci(ao) = cto ^ ^ ^^^'"'^ A is a constant independent of ao. Thus, near 00 the superpotcntial 
is given by 

W{u, ao) = ^ + Co + (ao - + a) u + ■ ■ ■ (20) 

00 u \ 2c-i J 

The result we have obtained depended crucially upon the assumption that c_i is not zero. If c_i = 0, the 
structure of the potential near 00 will be very different and the constraints given by Eqs. (|15p -(|17 p would 
no longer be valid. We would need to do the calculations again. This case is presented next. 

For c_i = 0, the superpotcntial near 00 is given by W{u,ao) — co + ciu + C2U^ + • • •. This leads to 
W'^{u,ao) = c§ + 2coCiM + {c\ + 2coC2)w^ + • • •, and = ciu^ + 2c2U^ + • • •• Since W'^{u,ao) and 

W'^{u,ai) must satisfy the shape invariance condition (|14p . matching the first two powers of u we get 

from the r° terms Co(ai) = Co(q!i) + R{ao) 

(21) 

from the r terms co(ao)ci(ao) = co{ai)ci{ai) . 

The first constraint of Eq. (f2T|) may be rewritten using ^ as c^{ao) + g{ao) ~ Co(ao + 1) +5(^0 + 1) which 
means that the quantity c^(ao) + 5(0^0) is independent of the argument aoi a-nd hence equal to a constant 
A. This leads to 



Co 



(ao) = ±V-5(«o)+A . (22) 



Now from the second constraint we see that the product co(ao)ci(ao) is also independent of the argument 
ao, and hence it must also be equal to a constant, which we denote by C,. Thus, from co(ao)ci(ao) = Cj we 
have 



ci (ao) = — ^ = ^ , f ^ ^ . (23) 
co(ao) ±V-5(ao) + A 



Thus, near 00, 



W{u,ao) 



± V-5(ao)+A ± 



v/-.g(ao) + A 



(24) 
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Thus, we have gathered significant knowledge about the structure of the superpotential near boundary points 
simply from the shape invariance. Recapitulating: near the origin, as given in Eq. (jl3p . the structure is 
1^ (r , ao) = - ^ + ^ + • • •. Near oo, as u ^ 0, there are two possible structures that depend on whether 
the superpotential has a singularity; i.e., whether the value of the coefficient c_i(q!o) is zero or non-zero. 
They are: 

r £^ + Co + (ao - S + a) u + • • • c_i(ao) + 
W{u,a^)^\ (25) 
±V-g("o)+A ± ^, " + ••• c_i(ao)=0. 

With the understanding about the structure of W{x, ao) at the end points of the domain, we now 
substitute these expansions of W into QHJ and obtain the implications on the QMFs, and hence on the 
eigenvalues of the hamiltonian. We first expand the momentum p(r) near the origin as: 



p{r, ao) = Pkr'' ^ + po f " + Pi + P2 



The derivative of the momentum p is then given by ^MLpl ^ + + Pi + 2p2r + ■ ■ ■ and the square 

of p is given by p^r, ao) = ^ + (^mI + (p2 + 2p_i p^) + ■ ■ ■. 

Collecting terms with various powers of r, we find that near the origin the combination p^ — p' is given 

by 

\P -^'^lorigin -2 + +(Po + 2p-lPl- Pi) + ■ ■ ' 

Substituting into the QHJ equation ([4]) near r = 0, we get 

p_i = -ao . (26) 

It is this term in the expansion of p that contributes to the contour integration aroimd the origin in the 
complex r plane. So, we do not need other coefhcients. 

To determine the leading order behavior of the momentum function at cxd, let us first consider the case 
that ci{ao) 7^ 0. We expand the momentum function p (i , ao) = p(it, ao) in powers of u: 



p{u,ao) = y^VkU = It pou° + Piu + P2U +■■■ 

fc=-i 

This leads to the following expression for the combination p^ — p': 

{P'-P')L= (p'+^'f) -p'iU-' + {2p-iPo)u-' + {p^ + 2p^ipi-p^i)+--- (27) 

Substituting into Eq. (g]), + = W^^ + - S , we obtain 

(P-i)^ , 2poP~i , /~2 , o~ ~ ~ ^ , ("^-i)^ , 2coc_i 

—^ + —— + {p, + 2p^ipi-p^i)+--- = _^ + __ (28) 

+ U + 2C., fao-'^)-c^i-E 



By matching the various powers of the variable u, we get p_i = c_i ; po = cq; and jjq + 2p_ipi — p_i 



2c_i ^cko — 2^~) ~ c_i — E. From the last equality, we get 2c_i ^ao 



g(Qo) 

2c_i 



= (2c_ipi), i.e., 



Pi — ao — ■ At 00, it is this last coefficient pi that contributes to the contour integration. 

We are now ready to determine the eigenvalues for the system using the information we have generated 
solely from shape invariance. The contribution to the contour integration is given by 



pdr 



Moving poles 



pdr 



r=0 



pdr 



pdr + (p pdu 
r=a Ju=0 ^ 



'{2iTip^i) + {2nipi) 



-27r i (— ao) + 2t: i 



ao 



9{ao)+E 



2 c- 



2TTi 



2ao 



giao) + E 



2 c- 



(29) 



Please note that due to the factor in the term §^^q[^ p) du, only the linear term of p contributes towards 
the integral and is given by (27ripi). 

It is worth noting that we cannot get the entire global behavior of the potential from expansions at the 
origin and 00. If the potential is a symmetric function of r, as is the case for the three dimensional harmonic 
oscillator, one gets an additional contribution [3, 5J. As one embeds the domain (0,r) in a two-dimensional 
complex plane, a mirror image of the potential function in the region TZeal(r) < also generates singularities 
for p(r) exactly at the same points as p(r) does in the region TZeal{r) > 0. As a consequence, we get 



pdr 



Moving poles 



pdr 



Tleal(r)>0 



pdr 



Tleal(r)<Q 



pdr 



r=0 



1 



pdu 



However, [§ pdr]^^^^^^^^^^ and [§ pdr]^^^^^^^^^ 
This leads to 



have identical value due to the symmetry of the potential. 



pdr 

From Eq. ([29]) . we then get 



Moving poles 



1 



2-Ki 



2ao- 



pdr 



g{ao) + En 
2c-i 



pdu 



= — 27rm 



-2t: in . 



(30) 



Solving for the energy, the Eq. (|29p gives En = 4c_i (n + ao) ~ 5('^o)- To determine the value of 17(00), we 
remember that for unbroken SUSY that the groundstate energy Eq — 0. Thus, substituting for n — 0, we 
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get g{a()) — 4c_iao, and hence g{an) = 4c_i {ag + n). Substituting for g{ao) and g{an) in i?„, we get 

i;„ = 4nc_i = g{an) - giao). (31) 

This is the desired result for any shape invariant problem. For the 3-dimensional harmonic oscillator, one 
would identify the coefhcient c_i with i u, and the energy would be given by 

En = 47ic_i = 2nu; , (32) 

as expected. However, the result of Eq. (j3ip is valid for any superpotential that is a) shape invariant and 
b) has singularities at both infinity and the origin. 

Now let us consider the case that the superpotential has no divergence at oo, i.e., c_i = 0. For this case, 
we will need to determine the structure of the quantum momentum function p{x) at oo (the structure at the 
origin remains the same as that for c_i 7^ 0.) Wc substitute the expansion of superpotential W{u,aQ) from 
Eq. (|25p into the Quantum Hamilton- Jacobi equation at 00, and we get 



+ J^2^ + (p2^2p_ipi-p_i) + --- = {-g{aa)+A-E) +2Cu + --- . 



Now equating the coefficients of various powers of u, we get p_i{ao) = 0; Po(Q!o) = ^5(^0) + A — ; 
Po{ao)pi{ao) = C; 6tc. From these, we obtain the following coefficients for p: 

c 



PQ 



= ±y/-g{aQ) + A - E; and pi = ± 



^-g{ao) + A-E 



The quantization condition yields: 
pdr 



Moving poles 



pdr 



r=0 



\ pdu 



M=0 



second = — (27rip_i) + {2'Kipi) 



= — 27ri (—ao) + 27r j 
= —2'Kin . 



Thus, we have — ao , ^ = n. Solving for En, 

±\/-g{ao)+^^-E 



±- 



C 



v/-5(ao) + A - £; 



En 



c 



+ g{ao) - A. 



(33) 



(34) 



{n + ao)2 

Since we are in the domain of unbroken SUSY, we must have E^ = 0. Hence substituting n = in Eq. (|34p . 
we get g{ao) = + A; i.e., 5(a„) = - (ao+n)'^ + ^- energy En is then given by 

1 1 



En = C 



= 9{oin) - giao). 



(35) 



_(ao)2 (n + ao)2_ 

Thus, we determined the energy of this system described by a potential that is singular at the origin and has 
no singularity at 00 (Coulomb- like) simply by employing the shape invariance condition. Once we identify 
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with e^/2 and ao with I + 1, this energy is exactly the same as that for the Coulomb potential. Once again, 
we state that Eq. (|35p is valid not only for the Coulomb potential, but for any shape invariant potential 
that has a fixed pole at the origin and none at the oo. 

Hyperbolic potentials: 

Now we consider superpotentials that arise as a linear combination of various powers of a basic exponential 
function e^. For example, the Morse, Eckart, hyperbolic Rosen-Morse potentials would fall under this 
category. We assume that our generic superpotential has the general form 

W{x, ao) = hj [e^'Y . (36) 

3 

This superpotential has singularities at ±oo. To analyze the nature of these singularities we first convert 
it into an algebraic form with a transformation r — e^. This change of the variable maps x G (— cxd,c») to 
r € (0, oo). We then subject the resulting superpotential to the shape invariance condition. 

Near negative oo; i.e., x — > — oo, and r — > 0, we assume the potential has the following structure: 

W(x,ao) = &-ir-^+6o + 6ir + ---. (37) 

Near positive oo; r — > oo, we express it in terms of a variable u = 1/r = e~^, where we assume it has the 
following structure: 

W{u,aQ)^c^iu-'^ +ca + ciu^ . (38) 

Now we will determine how shape invariance constrains various terms in these two expansions at the end 
points of the domain of r. 

As with the algebraic superpotential, the following procedure differs depending on whether 6_i and c_i 
are each zero or nonzero, making for four possible cases. We want to stress that we are considering only 
those shape invariant superpotentials that have at least one singularity at either the origin or oo. If the 
superpotential has no poles at these extremes, then it must have poles at "finite" fixed-points. These are 
potentials of trigonometric type, and as noted earlier, will not be considered in this paper. 

It is not hard to show that the case 6_i — 0, c_i 7^ is identical to the case 6_i 7^ 0, c_i = 0. Hence, 
we need consider only two cases: &_i 7^ 0, c_i = and &_i 7^ 0, c_i 7^ 0. However, it turns out that 
for (fe_i 7^ 0, c_i 7^ 0), the quantization condition § pdx ~ — 27rm, does not yield an equation involving 
energy, and hence energy cannot be determined from singularities at boundary points alone. The QMF in 
such a case must have poles at points other than the origin and 00 and, as stated earlier, such cases will 
not be considered here. It is important to note that there are no exponential type known shape invariant 
potentials in quantum mechanics for which the superpotential has a singularities for both x — ±00. Thus, 
in the following, we will only consider the case where 6_i 7^ and c_i =0. 
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Near r = 0, the superpotential is given by 

W{r, ao) = ^^i^ + 6o(ao) + h{ao)r + ■■■ (39) 
r 

From W, we get ^ = =r^^ + hr, and = '-^ + + fcg + 2b^ih + ■■■■ Hence, due 

to the above change of variable, the shape invariance condition now takes the form 

WHao) + = W^a,) - + i?(ao). (40) 

ar ar 

Substituting the expansion of W into Eq. (I40|) . and matching various powers of r, we get the foUowing set 
of difference equations: 

from the r^^ terms 6^]^(ao) = b^i{ai); 

from the r^^ terms 26_i(ao)&o(«o) — &-i(ao) = 25_i(Q;i)5o(ai) + 6_i(q;i); (41) 

from the r° terms ^o('^o) + 26_i(ao)6i(Q;o) — feo('^i) + 26-i(ai)6i(Q!i) + i?(ao). 

From Eq. (|4ip . the first constraint implies that 6_i(q;o) does not depend upon its argument; i.e., it is a 
constant. We denote this constant by The second constraint we get from matching the coefficients of 
gives 6o(Q!i) = ^o(<^o) ~ 1- Substituting ai = ag + 1, we get the difference equation 

bo{ao + 1) = 6o(ao) - 1, 

whose solution is 6o(q;o) = — ao + C, where C is a constant. Now, from the last constraint of the set: 6Q(ao) + 
26_i6i(ao) — ^0(0^1) + 26_i6i(q;i) + R{ao), we have 

(-ao + Cf + 26_i6i(ao) = (-ai + C)^ + 26_i6i(ai) + g{ai) - g{ao). 

This equation can be written as 

, , g(ao) + (-gp + C)' , , g(ai) + (-Qi + C)' 
61 (ao) + = 61 (ai) + , 

which shows that the expression 61 (ao) + 9(aa)+{-ao+C) independent of its argument ao, hence we set it 
equal to a constant D. Thus, we get 

We carry out a very similar analysis in the region r 00, where the superpotential has the form given in 
Eq. ([55]) : i.e., W{u, ao) = co(ao) + ci(ao) u + ci(ao) • • •, and we get 

(co(ao))^ = (co(ai))^ + i?(ao) 

Writing this equation as (co(ao))^ + g{ao) — (co(ai))^ + g(ai), we see that the expression (co(ao))^ + g{ao) 
is independent of ao. We set it equal to a constant 8. Thus, we get 

co(ao) = ±V6 - g{ao)- 
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Thus, to summarize, the structure of W at two end points are given by 



r 



D- 



2b-i 



2 \ 1 



and 



W{u,ao) = (±\/0 - .g(ao)) + ci{ao) u + ci{ao) ■ 



We will now substitute these two forms of the superpotential into the QHJ equation to determine the 
analytic structure of p and from it the eigenvalues of the shape invariant system. 
Near r — 0, the QHJ equation is given by 



2 w2/ ^ dW{r,ao) 

p —r— = W (r, aoj — r 
dr 



dr 



-E, 



where W{r, ao) is given by W{r, ag) = — + (— ao + C) 



D - 



) + (-ao+C)- 



26- 



quantum momentum hmction as p{r) = p-ir ^ + Po + Pif + ' ' ') and substituting it into Eq. 

p\ = b\; 2p_ipo - P-i = 2fe_i6o(ao) - etc. 



(43) 

Expanding the 
we get 

(44) 



From the first equality in above equation, we get p_i = ±&_i. However, in the limit £' — > 0, we must 
have p —^ W. This implies p-i — Substituting this value of p_i in the second equality, we get 

Po — ^o(q^o) = ^c^o + C. As we will soon see, we need go no further, pq is all we need to determine the 
contour integral around r = 0. 

Near u — 0, i.e., r oo, we expand the QMF as p{r) — p^iu~^ +Po + Piu + ■ ■ ■. We substitute this 
expansion for p{r) and the superpotential W{u,ao) = ^±^^6 — g{ao)j + ci{aa) u + ci{ao) u- 
QHJ equation 



into the 



p^+u^ = WHao)+u^^^-E. 
du du 



(45) 



Now comparing various powers of u, we get p_i — 0, (po)^ = ^±^0 — ^(ao)^ — E, etc. Thus, we get 



Po = ±^0 - 5(ao) - E. 

The quantization condition in this case is 

p dx - 

x — ~oo 



p dx = 2TTin ; i.e., 



r=0 



(?) 



dr 



u=0 



Po 



du = 271111 



(C - ao) ± a/6 - g(ao) - En = n , 



which yields. 



-B„ = 6 - (C - ao - n) - g{ao) 
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Since we assume that SUSY is unbroken, we must have Eq = 0. This imphes g{ao) = — (C — ao)^, and 
hence, g{an) = O — (C — (cto + n))"^ = Q — {C — aQ — Thus, we have 

En = {C - aof - (C - ao - nf = - g{ao) . (46) 

If we identify C — ag with the parameter A of Morse superpotential W = A — Be^^, we get the energy 
En = A? — {A — n)^, which is exactly what is expected for the Morse potential. 

Conclusion: One of the current authors (AG) and his collaborators had previously connected Quantum 
Hamilton- Jacobi Theory with supersymmetric quantum mechanics, and had shown that the quantum mo- 
menta for a shape invariant system can be derived recursively, analogous to SUSYQM. In this paper, we 
have further explored the impact of shape invariance on Quantum Hamilton- Jacobi theory. In particular, we 
showed that the shape invariance condition provides sufficient information about the singularity structure 
of QMF's to determine the eigenvalues of the system. 

There are several directions in which this work could be extended. We only considered translational 
shape invariant potentials; that is, those for which partner potentials are connected by a shift in parameter 
(fli-i-i = fli + 1). This leaves out the very important case where parameters in partner potentials are 
multiplicatively connected, i.e., a^+i = gci, where q e (0, 1). In addition, we only considered the potentials 
for which all fixed poles were at the origin or at oo or both. Potentials with finite domain, and those for 
which the superpotential has poles at points other than or/and co can also be explored. Recently, the 
work in SUSYQM has been extended to fractional supersymmetric quantum mechanics [121 [E] . One of the 
benefits of this formalism is that several solvable potentials are generated from one generalized potential for 
varying values of a parameter. The solvability in this case is described by a generalized shape invariance 
condition that is somewhat more involved than the additive shape invariance considered here. It would be 
worth exploring if the present formalism can be extended to fractional SUSYQM. 

We would like to thank Prof. John Dykla, Prof. Constantine Rasinariu and Walter Moore for discussions. 
We would also like to thank Prof. Jcffry Mallow for a careful reading of the manuscripts and providing many 
helpful suggestions. We would also like to thank one of the referees for several helpful suggestions. 

References 

[1] Cooper F, Khare A and Sukhatme U, Phys. Rep. 251 267, 1995; Cooper P., Khare A. and Sukhatme 
U., 2001 Swper symmetry in Quantum Mechanics (Singapore: World Scientific); G. Junker, 1996, Super- 
symmetric Methods in Quantum and Statistical Physics, Springer, Berlin. 

[2] Dutt R., Khare A. and Sukhatme U., Am. J. Phys. 56 163, 1988. 

[3] Leacock R. A. and Padgett M. J., Phys. Rev. Lett. 50 3, 1983; Phys. Rev. D 28 2491, 1983; Phys. Rev. 
A 33 2775, 1986; Am. J. Phys. 55, 261, 1986. 



13 



[4] Gozzi E., Phys. Lett. B 165, 351, 1985. 

[5] Bhalla R. S., Kapoor A. K. and Panigrahi P. K., Am. J. Phys. 65 1187, 1997; Mod. Phys. Lett. A 12 
295, 1997; Geojo K G, Sree Ranjani S. and Kapoor A. K., J. Phys. A: Math. Gen. 36, 4591-4598, 2003. 

[6] Rasinariu C, Dykla J., Gangopadhyaya A., Mallow J., Phys. Lett. A 338 197, 2005. 

[7] S. Sree Ranjani, Geojo K.G., Kapoor A.K., Panigrahi P.K., Mod. Phys. Lett. A 19, 1457-1468, 2004. 

[8] S. Sree Ranjani, Kapoor A.K., Panigrahi P.K., Int. Jour, of Theo. Phys., 44, 1167-1176, 2005; Ann. of 
Phys., 320 164-174, 2005. 

[9] A. Gangopadhyaya and U. P. Sukhatme, Phys. Lett. A 224, 5 (1996); U. P. Sukhatme, C. Rasinariu, 
and A. Khare, Phys. Lett. A 234, 401 (1997); Sukhatme UP, Rasinariu C, Khare A. and Gangopad- 
hyaya A, Supersymmetry and Integrable models, Springer, LNP 502, page 369-379. Proc. of a workshop 
at University of Illinois at Chicago, June 1997. 

[10] A. Gangopadhyaya and U. P. Sukhatme, Phys. Lett. A 224, 5 (1996). 

[11] Gangopadhyaya A., Panigrahi P. K. Sukhatme U. P., Jour, of Phys. A-Mathematical & General. 
27(12):4295-4300, (1994). 

[12] M. Daoud and M. Kibler, Phys. Lett. A 321, 147 (2004). 

[13] M. Daoud and M. Kibler, Jour. Math. Phys. bf 47, 122108 (2006). 



14 



